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^ Abstract 

Bilayer graphene in a magnetic field supports eight zero-energy Landau levels, which, as a tunable band gap develops, evolve into 
two nearly-degenerate quartets separated by the band gap. A close look is made into the properties of such an isolated quartet of 
^ pseudo-zero-mode levels at half filling in the presence of an in-plane electric field and the Coulomb interaction. The pseudo-zero- 
/— s mode levels turn out to support, via orbital mixing, charge carriers with induced electric dipole moment, which lead to characteristic 
l— ~ ' collective excitations, pseudospin waves, with some controllable features; possible experimental signals are discussed. 

t-H Key words: Graphene, bilayer, quantum Hall effect, collective excitations 
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1. Introduction 

Graphene attracts a great deal of attention, both experimen- 
tally and theoretically yj, |2J], for its exotic electronic trans- 
port. Of particular interest is bilayer graphene fl, which has a 
unique property that the (conduction-valence) band gap is con- 
trollable |4] by use of external gates or chemical doping. 

A particle-hole symmetric tower of Landau levels and the 
presence of some zero-energy levels are features specific to 
graphene in a magnetic field. Bilayer graphene supports eight 
such zero-energy levels, which, as a tunable band gap devel- 
ops, evolve into two nearly-degenerate quartets separated by 
i the band gap. 

In this paper we would like to report some unusual proper- 
ties of such an isolated quartet of pseudo-zero-mode levels, es- 
pecially, coherence and collective excitations in the presence of 
an in-plane electric field Ey and the Coulomb interaction 
We focus on revealing further controllable features in bilayer 
graphene and point out the following: 

(1) The pseudo-zero-mode levels, especially at half filling, 
support, via orbital level mixing, charge carriers with a nonzero 
electric dipole moment. As a result, fine splitting of the pseudo- 
zero-mode quartet is also externally controlled by an in-plane 
electric field or by an injected Hall current. 

(2) The interplay of the Coulomb interaction and an exter- 
nal field leads to rich spectra of collective excitations, (orbital) 
pseudospin waves. 

i 

2. Pseudo-zero-mode levels and dynamics 

Monolayer graphene supports as charge carriers massless 
Dirac fermions with a linear dispersion. In bilayer graphene 
interlayer coupling modifies the intralayer "relativistic" spectra 
to yield, in the low-energy branches, massive quasiparticlesH, 
which, in a magnetic field B, lead to a tower of Landau levels 
\n,y ) of energy 

e„ = sgn[«]<j*VM(M-l)> (1) 
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Figure 1 : Pseudo-zero-mode Landau levels 



labeled by integers n = +0, ±1, ±2, . . ., the center coordinate 
y = ( 2 p x with the magnetic length I = 1/ yfeB and the charac- 
teristic cyclotron energy u)^ 1 ~ 4 x B[T] meV (with B in tesla). 

The \n\ = and \n\ = 1 levels have zero energy, and, with a 
band gap U > 0, evolve into two sectors of nearly-degenerate 
pseudo-zero-mode levels (separated by a band gap and in val- 
ley), with spectra 



f«=±o 



±±t/, e„ =± i =±($U- 



2m). 



(2) 



where m w 0.006 x B[T] U <sc U denotes a fine splitting within 
each sector. 

The presence of the zero-energy levels and their double-fold 
degeneracy (per spin and valley) have a topological origin. This 
degeneracy is affected by a band gap U and an external field Aq. 
Accordingly, the pseudo-zero-mode levels have an intrinsic ten- 
dency to be degenerate, but its fine structure depends sensitively 
on the environment. 

The purpose of the present paper is to study such control- 
lable features of the pseudo-zero-mode quartet in the presence 
of external fields and Coulomb interactions. We are particularly 
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interested in the properties of this quartet at half filling, where 
mixing of the pseudo-zero-modes may take place. For defi- 
niteness we focus on the n = (0 + , 1) sector at one valley, i.e., 
around filling factor v = 2 (or v = 1 when the spin is resolved), 
and ignore the presence of other levels which are separated by 
relatively large gaps. 

The pseudo-zero-mode levels of our interest contain two dif- 
ferent Landau orbitals n = 0+ and n = 1 (per spin). Let us 
put these modes, ^o + (yo, f) and if/\(yo, f), into a two-component 
spinor *F = (^o + ,^i) t , and define the pseudospin operators in 
the (0 + , 1) orbital space, 



(3) 



with Pauli matrices <j° and <x° = 1; y p — e~ f2p2 / 4 . The charge 
density p_ p = J d 2 xe' px t/rty projected to the (0+, 1) sector is 
thereby written as 



p_ p = 2(w° p S°_ p + w 3 p S 3 _ p + w p Si p + w 2 p S 2 _ p ), 



(4) 



where w2 = 1 - £ 2 p 2 /4, wl = £ 2 p 2 /4, wi = iip y /^2, and 



w\p) = ilp x /^2. 



The relevant part of the bilayer Hamiltonian governing the 
electrons in the (0+, 1) sector is written as 



AH = - J d 2 x eA p + jdyom >PV 



(5) 



where Aq stands for an external potential and m a 0.006 x 
B[T] U for the fine level splitting. One can equally write 



AH = l^f^S^, 



(6) 



el 



(P°,P 3 ) * (-m-M ,m),(P 2 ,P 1 ) = -(£ I ,£ J ), (7) 

V2 

where Ey = (E x ,E y ) = -8 x Aq denotes the in-plane electric 
field. For the Coulomb interaction one can retain only the in- 
tralayer interaction since the pseudo-zero-modes at each valley 
essentially lie on the same layer, 



PpPp 



(8) 



where v p = 27ra/(eblpl) with a = e 2 /(4-neo) « 1/137 and the 
substrate dielectric constant e^. 



3. Pseudospin textures 

In this section we study the properties of the pseudo-zero- 
mode levels at half filling, using the projected Hamiltonian H = 
AH + H c , with focus on orbital mixing of the zero-modes. 

Let us suppose that the half-filled state |G) is given by a con- 
figuration where the pseudospin points in a fixed direction in 
pseudospin space, i.e., (G\S^ =0 \G) = jN e n a and n a n a = 1 
with the total number of electrons N e = 2{G\S° =0 \G). For n a 
we use the parametrization n 1 = sin# cos</>, n 2 = sin6> sin0, 



Note that « 3 = 1 corresponds to the filled n = 0+ level while 
n 3 = -1 represents the filled n = 1 level. The direction n = 
(n l ,n 2 ,n 3 ) would, in general, vary in response to the external 
field Aq, and, as n tilts from n 3 = ±1, the n = 0+ and n 
I levels start to mix. For selfconsistency we assume that Aq 
represents a uniform in-plane electric field Ey = -8 x Aq and that 
it leads to a homogeneous state \G) of uniform density p = 
v/(2n{ 2 ) (with v = 2 for the spin-degenerate v = 2 state). We 
thus consider all configurations with (G\S^ =0 \G) = \ N e n a and 
single out the ground state or the associated n a by minimizing 
the energy (G\H\G). 

The calculation of the Coulomb energy (G\H C \G), in par- 
ticular, requires the knowledge of pseudospin structure factors 
(G\SpSq\G), which, for the present half-filled state with pseu- 
dospin oc n a , are given in Ref. [5]. The energy is conveniently 
written in x space as {G\AH + H C \G) = po J d 2 ^.^^ with 



74ir = -M + E(6) - iVi, 

E(0) = m (cos - 1) + 6 sin 6 + ^ V x (1 - cos 6») 2 , (9) 



where £ = (ell V2) Ey • fty and fly = (h x , h y ) = (sin 0, cos <p) is 
an in-plane unit vector; 

(10) 



V x = yjn/2a/(e b l) 

sets the scale of the Coulomb exchange energy. In the present 
notation ?4ff stands for energy per electron in state |G). 

Note first that the half-filled pseudo-zero-mode state has an 
in-plane electric dipole moment 



el 

vi 



(« ,n ) 



el 



sin 6 nil 



(11) 



of strength (e^/ V2) | sin 0\ per electron, proportional to the in- 
plane component (n 2 ,n l ) of the pseudospin. Mixing of the 
n = 0+ and n = 1 modes gives rise to this dipole. Note also that 
the Coulomb correlation energy oc V\ is symmetric about the n 3 
axis || B. It alone favors 6 = 0, the filled n = 0+ level, and varies 
continuously by an amount AE C = (1/8) V] as n sweeps in 
pseudospin space. The Coulomb interaction thus significantly 
enhances the splitting of the pseudo-zero-mode levels. 

Obviously the energy is lower when when ny || Ey, or 
S = el\E\\\/ y2. Accordingly it is convenient, without loss of 
generality, to suppose that the in-plane field Ey and ny lie along 
the y axis, & = elE y / a/2 > 0. With this choice the "1", "2" and 
"3" axes in pseudospin space coincide with the y, x and -z axes 
in real space, respectively. We henceforth adopt this choice and 
set (n',n 2 ,n 3 ) = (sin#, 0, cos#). 

Let us explore stable configurations of the half-filled zero- 
mode state in the next section and here study collective excita- 
tions over a given ground state |G)|„. We focus on a special class 
of low-energy collective excitations, pseudospin waves, that are 
rotations about the energy minimum |G)|„ in pseudospin space. 
Such a collective state is represented as a texture state 



|G> = e - ,0 |G>| n , 



(12) 



cos 9, with -n <0 <n and < <f> < n. 



where the operator O = £ p Tp 1 ^ p S a p locally tilts the pseu- 
dospin from n by small angle £2 p . The energy (G\H\G) = 



2 



{G\e i0 He-'°\G) of the excited state is a functional of Q p , and 
serves as an effective Hamiltonian. The associated Lagrangian 
is neatly written as 

po J d 2 x L eff = {G\e i0 (id, - H)e- i0 \G). (13) 

Actual calculation is made in a systematic way if one notes that 
the pseudospin operators S p in Eq. (0 obey the SU(2)xWoo al- 
gebra. The result is an effective Lagrangian of the form ^ 

U = ^(3,cD) 2 -i<D(M p ) 2 0, 

M p = 2 ^ + F p ){k] + G p )-\W p \\ 

k 2 = E"(6), j^ = -£/sin0, (14) 

where (F p , G p , W p ) ~ 0(V\) go to zero for p — > 0. The spec- 
trum M p of pseudospin waves is in general anisotropic in p and 
depends critically on the stable-state configuration n. In partic- 
ular, Mp^o = and Mp^ * Vi/2 + k 2 + k 2 . With 6» -> 
and 3 — > 0, our M p precisely reproduces an excitation spectrum 
derived in Ref. [6] by assuming spatial isotropy. 

Note here that the Coulomb interaction alone yields k 2 = 0. 
This implies that, unlike in ordinary bilayer quantum Hall sys- 
tems, there is no cost of interlayer capacitance energy. Coher- 
ence is thus easier to form in the pseudo-zero-mode sector of 
bilayer graphene. 

4. ground states and pseudospin waves 

Let us now examine orbital configurations of the half-filled 
pseudo-zero-mode state and the associated pseudospin waves. 
Note first that the Coulomb correlation energy 

AE C = i jf SL- » (2.2/e sc ) VbTTI meV (15) 

with a typical value ~ 4 sets the basic energy scale; here 
e sc effectively takes account of the effect of screening coming 
from vacuum (Dirac-sea) polarization [5], specific to graphene. 
In ratio the intrinsic zero-mode level gap 2m and the in-plane 
fieldfi = e£\E n | /a/2 read 

m/AE c ~ 3e sc x 1(T 3 yfB\J]U[meV], (16) 
8/AE c * 0.9 e sc x 10" 3 £[V/cm]/B[T]. (17) 

One has to look for the minimum of E(6), E m [ n = E(6 m i n ), 
to determine the orbital configuration of |G). Let us begin with 
the case where Ey is absent, (i) For m > AE C (though rather 
unrealistic), one finds E m i n at 6 — n, i.e., the filled n — 1 
level is realized. The pseudospin waves have a finite energy 
gap M p= o = 2 (m - AE C ) > and the spectrum is isotropic. 

On the other hand, (ii) for < m < AE C , one finds E^ at 
6 = +9 m i n \ G m i n varies from to n with increasing m. Here we 
encounter a somewhat strange situation: The rotational invari- 
ance (about the applied field B) is spontaneously broken and the 
pseudospin waves are gapless. 

This actually signals an instability. In this m range the as- 
sumed ground state |G), taken to be homogeneous in space, 
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Figure 2: Pseudospin-wave excitation gap M p= o, in units of meV and GHz, 
plotted as a function of the in-plane field |Bn| for some typical values of mag- 
netic field B[T] and band gap f/[meV]. Dashed curves refer to the cases where 
the screening effect is turned off, e sc — > 1 . 

acquires spontaneous in-plane electric polarization oc sinf3 m i n . 
Such electrically-polarized homogeneous configurations, un- 
less polarization is relatively weak, are unstable against local 
charge inhomogeneities and would decay into inhomogeneous 
configurations. We thus speculate that the half-filled state in 
the realistic < m < AE C range may form many domains about 
local charge excesses for stabilization. 

(iii) The in-plane field & = e€E y l \2 > tilts the pseudospin 
toward 6 = -n/2 and competes with AE C . The charge carriers 
thereby acquire a nonzero electric dipole moment oc sin 6 and 
the pseudospin waves always have a finite excitation gap. The 
potential instability of the texture state, mentioned in (ii), weak- 
ens and eventually disappears for large |E|||. In Fig. 2 we plot 
the excitation gap M p= o as a function of |E||| for some typical 
values of B[T] and m oc {/[meV]; M p= o falls in the frequency 
range of microwaves. 

If E|| is sufficiently strong, a sizable gap M p= o would arise, 
leading to an incompressible v - 2 state. One could thereby 
observe the (spin-degenerate) v - 2 Hall plateau with a suitably 
strong injected current. Ey is practically strong for£ > 0.3 AE C , 
where M p= o = 2 « 2& + j AE C with k 2 m S + j AE C and 

It would be a challenge to directly detect the pseudospin- 
wave excitation gap by microwave absorption or reflection ex- 
perimentally. One may follow Ref. [5] and derive the optical 
conductance due to transitions within the (0+, 1) sector, 

Acr«(w) = (ve 2 IK)K 2 5(u)- M p=0 ), 

Ao- yy (oS) = {ye 2 IK) (cos O^f 6(oj - M p=0 ), (18) 

which, with disorder taken into account, is significantly peaked 
around u ~ M p= q. Note that the optical conductance is spatially 
anisotropic, with the larger component Acr xx (tS) perpendicular 
to E|| (or along an injected current). 
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5. Summary and discussion 



Zero-mode Landau levels, specific to graphene in a mag- 
netic field, deserve attention in their own right. In bilayer 
graphene they evolve, with a tunable band gap, into two quar- 
tets of nearly-degenerate pseudo-zero-mode levels, which, un- 
like in monolayer graphene, involve two different orbital in- 
dices n = 0, 1 . One would thus expect interesting quantum 
phenomena of orbital level mixing. 

In this paper we have studied the effects of an external field 
and the Coulomb interaction on the pseudo-zero-mode quartet. 
This quartet, especially at half filling, supports, via orbital mix- 
ing, quasiparticles with charge and electric dipole, which give 
rise to characteristic collective excitations, pseudospin waves. 

The pseudospin-wave excitation gap M p= o is generally small, 
reflecting the intrinsic degeneracy of the pseudo-zero-mode lev- 
els, but turns out to depend sensitively on an in-plane field Ey. 
This means that the gap M p= o is tunable by an in-plane field or 
by an injected current. 

An experimental signature of the field-induced gap is to 
observe the quantum Hall effect with an injected current; one 
would be able to resolve the v = ±2 Hall plateaus (or the spin- 
resolved v = +1 plateaus) using a suitably strong current. A 
direct study of the excitation gap M p= o and its field dependence 
by microwave absorption or reflection would also clarify the 
unique controllable features of the pseudo-zero-mode sector in 
bilayer graphene. 

This work was supported in part by a Grant-in- Aid for Scien- 
tific Research from the Ministry of Education, Science, Sports 
and Culture of Japan (Grant No. 21540265). 
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